Systematics of quark— ant iquark states: where are the lightest glueballs? 



V.V. Anisovich 

(Talk given at HADRON-2003, 31 August - 5 September 2003, Aschaffensburg, Germany) 

The analysis of the experimental data of Crystal Barrel Collaboration on the pp annihilation in 
flight with the production of mesons in the final state resulted in a discovery of a large number of 
mesons over the region 1900-2400 MeV, thus allowing us to systematize quark-antiquark states in 
the (n, M 2 ) and (J, M 2 ) planes, where n and J are radial quantum number and spin of the meson 
with the mass M. The data point to meson trajectories in these planes being approximately linear, 
with a universal slope. Basing on these data and results of the recent K-matrix analysis a nonet 
classification is performed. In the scalar-isoscalar sector, the broad resonance state /o(1200 — 1600) is 
superfluous for the qq classification, i.e. it is an exotic state. The ratios of coupling constants for the 
transitions fo — * tvtv, KK ,r]rj,r]rj' point to the gluonium nature of the broad state /o(1200 — 1600). 
The problem of the location of the lightest pseudoscalar glueball is also discussed. 

The search for exotic mesons should be based on the classification of (/(/-states. Exotic mesons are those which are 
superfluous for the qq systematics. The quark-antiquark systematics means: 

(i) classification of qq states as states located on the (n, M 2 ) and (J, M 2 ) trajectories, and 

(ii) determination of the quark-gluonium content of states from the analysis of the decay coupling constants, namely, 
hadronic and radiative decay couplings as well as weak ones. 

For the hadronic decay coupling constants, the most reliable information comes from the -ftT-matrix analysis. In 
addition, the if-matrix analysis allows us to study bare states (the states before the onset of the decay processes). 

1. Systematics of the gg-states on the (n,M 2 ) and (J, M 2 ) planes. 

The analysis of experimental data on the pp annihilation in flight with the production of mesons in the final state 
resulted in a discovery of the large number of mesons over the region 1900-2400 MeV [1] . This allowed us to systematize 
quark-antiquark states on the (n, M 2 ) and (J, M 2 ) planes. The data point to almost the linear meson trajectories on 
these planes, with a universal slope [2]. 

In Fig. 1, one can see the (n,M 2 ) trajectories for the (I = 1) states, which are drown for the a\- and a3-mesons 
(Fig. la), 7T-, 7T2- and 7r4-mesons (Fig. lb), b%- and 63-mesons (Fig. lc). All these trajectories reveal linear behaviour, 
such as 

M 2 = M 2 + /i 2 (n-l), 

with fj 2 ~ 1.2 GeV 2 ; Mq is the mass of the ground (basic) state, n = 1. The pion, being beyond the trajectory, is 
an exception, that is not a surprise, for the pion is a special particle in certain respect. In the classification, all these 
mesons should be treated as gg-states. Using the spectroscopy notations for gg-states, 2S+1 Lj where S is the quark 
spin and L is their orbital momentum, we assign the trajectories to mesons as follows: 

ai (1230) trajectory: n 3 Pi gg-states, 03 (2030) trajectory: r£ 3 i<3gg-states; 

7r(140) trajectory: n 1 S , gg-states, 7r 2 (1670) trajectory: n 1 D 2 qq-st&tes, 7r 4 (2250) trajectory: ^(S^ gg-states; 
6i(1235) trajectory: n 1 Pi gg-states, 63(2020) trajectory: T^i^gg-states. 

In Fig. lc the state 61 (1640) is shown which was not discovered in the experiment but is predicted by trajecto- 
ries: the states we predict are denoted by open circles. 

The trajectories p, P3, ao, 02, a 4 demonstrate linear behaviour as well. The p-trajectories are shown in Fig. 2a: 

p(770) trajectory: dominantly n 3 Si qq{L — 0), p(1700) trajectory: dominantly n 3 Di qq(L = 4). 

The states 3 Si gg and 3 Di qq may mix with each other but considerable mass splitting of the p(770) and p(1700) 
states tells us that the mixing is not large. The and ao, &2, a 4 trajectories are shown in Figs. 2b and 2c: 

/33(1690) trajectory: dominantly n 3 Ds qq(L — 2), p3(2260) trajectory: dominantly n 3 G?, qq(L = 4), 

ao(980) trajectory: n 3 P Q qq(L — 1), 02(1320) trajectory: dominantly n 3 P 2 qq(L = 1), 

a2(2030) trajectory: dominantly n 3 F2 qq(L — 3), a 4 (2005) trajectory: dominantly n 3 Fn qq(L = 3). 



1 






2 



The ao(980)-resonance, which sometimes is discussed as a non-qq state, lays on the gg-trajectory, that is an argument 
in favour of its quark-antiquark origin. 

The figures 3 and 4 for the 1 = states confirm the linearity of trajectories on the (n, M 2 )-planes, with a universal 
slope /i 2 ~ 1.2 — 1.3 GeV 2 . In this sector, we face a doubling of trajectories due to the existence of two flavour states 
nn = (uu + dd)/y/2 and ss. 

In Fig. 3, one can see the /2(1285) and /2(1525) trajectories, which are dominantly the 3 P2 (nn + ss) states, and 
/2(2020) and /4(2020) ones, which are dominantly 3 F2 (nn + ss) and 3 F4 (nn + ss), correspondingly. Figure 3b 
demonstrates /o(980) and /o(1300) trajectories ( 3 Po (nn + ss) states), while figure 3c shows the trajectories w(780), 
0(1020) (dominantly 3 <Si (nn + ss)) and w(1640) (dominantly 3 Di (nn + ss)). 





FIG. 5. The it and ai trajectories on the (J,M 2 ) plane. 

Let me stress that / (980), / (1300) (the latter denoted as / (1370) in the compilation [4]), / (1500), / (1750), 
which are sometimes discussed as candidates to exotics, lay quite comfortably on the linear qq trajectories. The K- 
matrix analysis [5] gives us one more state in the (IJ PC = 00 ++ )-sector, that is, a broad resonance /o(1200 — 1600): 
just this state may be considered as exotics, a descendant of the lightest scalar glueball. It is discussed in the next 
section. The light a, if it exists, is also beyond the qq trajectories being in this way a candidate for exotics as well. 

In Figs. 4a,b,c, one can see trajectoties as follows: u 3 , (dominantely 3 D 3 qq), hi ( 1 P i qq), r] ( 1 S'o qq) and 772 
( 1 i?2 qq)- Let me emphasize that the situation in the ^ sector is more complicated than it is seen from Fig. 4c. On 
the one hand, the states ry(1440) and 7/(1295) lay good enough on the linear (n,M 2 ) trajectory. On the other hand, 
experimental indications for 7/(1295) are not convincing, and the resonance 7/(1440) reveals itself in different reactions 
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FIG. 7. The /a and trajectories on the (J, M 2 ) plane. 

with different masses: at 1390-1430 McV in the r]irir and ao(980)7r — * KKtt modes, while in the K*(890)K mode at 
1460-1500 MeV. Besides, in the reaction J/ip — * 777(1440), the resonance 7/(1440) is produced with a large background, 
that may tell us about the existence of a broad state at 1400-1700 MeV. Linear trajectories in the h sector predict 
two qq states in the neighbouring region, namely, ry(1700) and 7/(1820), but these resonance were not seen till now, 
because, with available data this mass region is difficult for the study. 

Now let us discuss the location of mesons on the (J, M 2 )-planes. The pion and a% trajectories are shown in Fig. 
5a, b together with their daughter ones. The trajectories for different parities are degenerate — this fact is illustrated 
by Fig. 5c, where the combined presentation of 7r and d\ trajectories is given. 

The (J, M 2 ) trajectories for p and ai depicted in Fig. 6 provide us unambigous information on the ao(980) 
resonance. One can clearly see that ao(980) lays on the linear daughter trajectory — this is obvious from their 
combined presentation in Fig. 6c. Supposing that ao(980) is non-gg meson, one should expect another ciq state, 
dominantly qq, around 1 GeV. However, additional ao state is definitely excluded by the experimental data. 

The (J, M 2 ) trajectories for the f% and u> resonances as well as for their daughter ones are shown in Fig. 7, and the 
combined presentation of Fig. 7c demonstrates the Jo (980) resonance laying on the daughter qq trajectory. 

At last, figure 8 shows us the trajectories for the iVmeson sector, the kaons with positive and negative parities lay 
on the degenerate trajectories. The k meson discussed as a plausible + state with the mass ~ 900 MeV does not 
belong to linear trajectory, so it should be considered as an exotic state. 



4 



Table 1. Nonet assignment of the qq states, n — 1 and 2. 



qq 

states 


n = 1 


n = 2 


1=1 1=0 1=0 1=1/2 


1=1 1=0 1=0 1=1/2 


<->ol.u ) 
3 5i(l-) 


p(770) w(780) 0(1020) tf*(890) 


7T^loUUJ ?7^1Zyoj 7/^144UJ i\^14DUj 

(0(1450) w(1420) 0(1657) #i(1680) 


3 Po(0++) 
3 Pi(l++) 

3 P 9 (2++) 


a (980) / (980) / (1300) X (1430) 
ai(1230) /i(1285) /i(1510) #i(1400) 

DiyLZoo) fli{±±i\J) UiyLovU) i\iyLZi\J) 

a 2 (1320) / 2 (1285) / 2 (1525) X 2 (1430) 


a (1520) / (1500) / (1750) tf (1850) 
ai(1640) if 1 (1780)? 

0i^1D4UJ: fti^lOUUJ: fliyiiou): i\iy±Vo\J) 

a 2 (1660) / 2 (1640) / 2 (1790) ^(1980) 


S D X {1— ) 
^2(2—) 

3 ^ 3 (3-) 


/3(1700) w(1640) ^(1680) 

X 2 (1800) 

*w„(~\(\7ft\ n^t'i&Atw n„(~\%(\r\\ je n h ^sn"! 

7T 2 ^1U/Uy // 2 ^lu^tOj // 2 ^±OUUy J\ 2 ^10OU7 

p 3 (1690) w 3 (1670) 3 (185O) # 3 (1780) 


(0(1990) w(1920) 

#2(2170)? 

7T 2 ^ZUUOy // 2 ^ZUOU^ /\ 2 ^ZUOuj: 

p 3 (1980) w 3 (1960) 3 (215O)? if 3 (1780) 


3 F 2 (2++) 

^3(3++) 

^3(3+-) 
3 F 4 (4++) 


a 2 (2030) / 2 (2020) 

a 3 (2030) i*f 3 (2030)? 
63(2020) i^ 3 (1960)? 
a 4 (2005) / 4 (2020) ^(2045) 


a 2 (2310)? / 2 (2290) 

a 3 (2275) X 3 (2320) 
63(2245) if 3 (2220)? 
a 4 (2255) / 4 (2300) 



Table 2. Nonet assignment of qq states, n — 3 and 4. 



qq 

states 


71 = 3 


n = 4 


1=1 1=0 1=0 1=1/2 


1=1 1=0 1=0 1=1/2 


L S (0- + ) 
3 5i(l~) 


tt(1800) 77(1700)? 7/(1820)? #(1830) 
p(1830)? w(1800)? 0(1950)? 


tt(2070) 77(2010) 
p(2150) w(2150) 


a Po(0 ++ ) 
3 Pi(l ++ ) 
1 P 1 (1+-) 
3 P,(2++) 


a (1830)? / (2005) / (2105) 
ai(1960) 

61 (1970) /ii(2000) /ii(2120)? 

a 2 (1950) / 2 (1950) / 2 (2110)? # 2 (2370)? 


a (2120)? / O (2230)? / (2330) 
ai(2270) 

6i(2210) /ii(2270) 

a 2 (2255) / 2 (2210) / 2 (2370)? 


8 A(1~) 
3 D 2 (2— ) 

1 D 2 (2-+) 
3 ^ 3 (3-) 


p(2285) w(2295) 

#2 (2400)? 

tt 2 (2245) 77 2 (2250) # 2 (2250) 
p 3 (2300) w 3 (2300)? 





Summing up the results of meson systematization on the (n, M ) and (J, M 2 ) trajectories, one can state that the 
resonances / (980), / (1300) (/ (1370) in [4]),/ (1500) and / (1750) (/o(1710) in [4]) are located on quark-antiquark 
trajectories. They originate from the standard quark-antiquark states. 

Now we can assign mesons to the flavour multiplets (see Tables 1, 2 for the multiplets with n = 1,2,3,4). The 
interrogation sign marks the states, which were not discovered by the experiment but are predicted by linear (n, M 2 ) 
and (J, M 2 ) trajectories. The empty places are left for the states, which were not discovered and cannot be predicted. 
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FIG. 8. Trajectories on the (J, M 2 ) plane, Tf-meson sector. 
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One may conclude that there are extra 00 ++ states with respect to the qq systematics on the (n,M 2 ) and (J, M 2 ) 
trajectories: 

1) the broad state /o(1200 — 1600) (it follows from the #-matrix analysis that it is the descendant of the glueball), 

2) the light a-meson, / (300— 500) (if it exists), 

3) K-meson, # (700 - 900) (if it exists). 

In the pseudoscalar-isoscalar sector (IJ PC = 00 h ), the situation in the mass region 1400-1800 MeV is rather 
uncertain, and one cannot exclude the states existing here, which are superfluous for the qq trajectories. 



2. The #-matrix analysis of the (IJ PC = 00 ++ )-wave. 

The trajectory assignment does not specify the content of the state. Each meson/resonance is a mixture of different 
components, and their wave functions are the Fock columns. For example, a and / may be considered as: 

qq 

qqqq 

7T7T, ##, rjrj, . 



ao 



qq 

qqqq 

777T, ##,... 



, fo 



99 



In bare states, determined by the #-matrix analysis, the long-range hadronic components are excluded that makes 
possible more reliable determination of the quark-antiquark and gluonium (gg) components. 
In the #-matrix analysis the partial wave amplitude reads: 

A(s) = K(s) \l - ipK(s) 

For the 00 ++ wave, which was analysed in [5], the following five channels were taken into account: 1 = irir, 2 = K #, 
3 = r]T], 4 = r/i] 1 , 5 = 7T7T7r7r, with the two-particle phase spaces determined as 



1 ^ 


#11, #12,- 




Pl,0, ... 


, # = 


#21, #22, ••• 


, P = 


,(>2,- 



Ami 



4m 2 . 



Pkk = 

v a 

while Ptttttt 7t was considcrd as cither two-cr or two-p phase space factor. 

The fitting parameters are the #-matrix elements, which are represented as the sums of pole terms, g^ 9^ /{p 2 L — s), 
and a smooth s-dependent term f a b(s): 



K ab 



(n) (n 



pI- s 



M n being the masses of bare states and g a , g b their couplings. 

The combined #- matrix analyses of the spectra have been performed for the mass interval 280 < M < 1950 MeV 
by including the following final states: 

7T7T, 7777, ##, 7777', 7T7r7r7r (pp, 00) [5], 

TFT], ## [6], 

#tt; [7]. 



1 = 
I = 1 
/= i 



Note that these combined #- matrix analyses have their predecessors published in [8]. The necessity of a combined 
analysis owes to the existence of large interference effects "resonance-background" as well as the effects associated 
with the resonance overlapping. In a situation of such a type, only a combined fitting to a large number of reactions 
allows one to expect reliable results. 

Previous analysis [6] carried out in 1997-1998 was based on the experimental data as follows: 

(1) GAMS data on the S'-wave two-meson production in the reactions irp — ► 7r°7r°n, -q-qn and rjij'n at small nucleon 
momenta transferred, |f| < 0.2 (GeV/c) 2 [9,10]; 

(2) GAMS data on the tttt S'-wave production in the reaction np — > 7r°7r°n at large momentum transfers squared, 
0.30 < |i| < 1.0 (GcV/c) 2 [9]; 

(3) BNL data on the reaction ir~p — > KKn [11]; 

(4) CERN-Miinich data on tt+tt" tt+tt" [12]; 

(5) Crystal Barrel data on pp (at rest, from liquid #2)— » 7r 7r°7r°, 7r°7r°?7, tt°i]T] [13,14]. 

Now the experimental basis has much broadened, and additional samples of data are included into the analysis [5] 
of the 00 ++ wave as follows: 
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(6) Crystal Barrel data on proton-antiproton annihilation in gas: pp (at rest, from gaseous H 2 ) — > 7r°7r 7r°, 7r°7r°?7 
[15]. 

One should keep in mind that in liquid hydrogen the pp annihilation is going dominantly from the S'-wave state, 
while in gas there is a considerable admixture of the P-wave, thus giving us an opportunity to analyse the three-meson 
Dalitz plots in more detail. 

(7) Crystal Barrel data on proton-antiproton annihilation in liquid: pp (at rest, from liquid H2)— > 7r + 7r~7r°, K + K~tt°, 
K s K s n°, K+K s tt- [15]; 

(8) Crystal Barrel data on neutron-antiproton annihilation in liquid deuterium: np(at rest, from liquid D 2 )^> ir°n°n~ , 
tt-tt-'tt+, K s K-tt°, K s K s n- [15]. 

These data allowed us to perform more confident study of the two-kaon channels as compared to what had been 
done before. This is important for the conclusion about the quark-gluon content of scalar-isoscalar / -mesons under 
investigation. 

(9) E852 Collaboration data on the irir S'-wave production in the reaction ir~p — > 7r°7r°n at the nucleon momentum 
transfers squared < |t| < 1.5 (GeV/c) 2 [16]. 

Experimental data of the E852 Collaboration on the reaction n~p — > 7r°7r°n at pi a b = 18 GcV/c [16] together with 
the GAMS data on the reaction n~p — > 7r°7r°ri at pi a t = 38 GeV/c [9] give us a solid ground for the study of the 
resonances /o(980) and /o(1300), for at large momenta transferred to the nucleon, |i| ~ (0.5 — 1.5) (GeV/c) 2 , the 
production of resonances is accompanied by a small background, thus allowing us to fix reliably their masses and 
widths. 

The most important ingredients of the new analysis [5] are: 

(i) the study of K ^-spectra that led to the determination of the flavour-octet and flavour-singlet components; 

(ii) the study of / (1300) produced without background: irN — > ttttN at large 

Figure 9 taken from [5] demonstrates the complex-M plane for the 00 ++ sector. Here the masses and total widths 
of resonances are determined by the position of amplitude poles, det\l — ipK\ — 0, the decay couplings are determined 
by the pole residues. 

The movement of poles in the complex-M plane for the states / (980), / (1300), / (1500), / (1750), / (1200- 1600) 
with a uniform onset of the decay channels is shown in Fig. 10. Technically, to switch on/off the decay channels for 
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FIG. 9. Complex-M plane in the (IJ PC = 00 ++ ) sector: masses and total widths of resonances are determined by the 
position of the amplitude poles det|l — ipK\ = 0, decay couplings are determined by the pole residues. Dashed line encircles 
the part of the plane where the if-matrix analysis reconstructs the analytical /^-matrix amplitude: in this area the poles 
corresponding to resonances /o(980), /o(1300), /o(1500), /o(1750) and the broad state /o(1200 — 1600) are located. On the 
border of this area the light cr-meson denoted as / (450) is shown (however, we have also rather good solution without a- meson 
pole). 

the K -matrix amplitude one should substitute in the if- matrix elements — > £n( x )3a^ and / a & — > £f(x)f a b, where 
the parameter-functions £ n (x) and £/(x) satisfy the following constraints: ^ n (0) = £f{0) = and ^ n (l) = £/(l) = 1, 
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and x varies in the interval < x < 1. Then, at x = 0, the amplitude A turns into the if -matrix, A(x — > 0) — > If, 
and the amplitude poles occur on the real axis, that corresponds to stable /Q° re -states. 
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FIG. 10. Complex M plane: trajectories of poles corresponding to the states /o(980), /o(1300), /o(1500), /o(1750), 
/o(1200 — 1600) within a uniform onset of the decay channels. 



In Fig. 10, one can see gradual transformation of bare states into real mesons as follows: 
/ bare (700 ± 100) -» / (980), / o jare (1220±40) -» / (1300), 



/ ( 6ore (1230±40) 
.f bare (1800±40) 



/o(1500) 
/o(1750) 



/ 6are (1580±40) -» / (1200 - 1600) 



Although in the K-matrix analysis the pole position of the broad state /o(1200 — 1600) is defined with large errors, 
this pole is necessary in all the if-matrix solutions. Also in all solutions the requirement of the factorization of 
coupling constants is fulfilled: 
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The onset of the decay channels can be illustrated with an example of the / -levels in the potential well (Fig. 11): 
stable levels correspond to bare states (Fig. 11a), while overlapping resonances correspond to real mesons (Fig. lib). 

Figure 12a demonstrates that bare states also form linear trajectories. Here one can see two /g are trajectories as 
well as Og are and K^ re ones: all these trajectories have almost the same slopes as the trajectories of real resonances, 
which are also shown in Fig. 12b for the comparison. 

To fix the nonet of bare states by using hadronic decay processes two parameters are needed only, namely, g, ip, 
where g is a universal decay coupling and <p the mixing angle for nh and ss components in /o° re (l) and /Q° re (2). 
Decay couplings depend also on the suppression parameter for the strange quark production probability A — 0.6 ±0.2, 
but its value is fixed by other reactions, for example, see [17]. 

Let us emphasize that these two parameters, g and (p, allowed us to describe ten decay reactions such as 

/o(l) -> it-it, KK,rrn, / (2) -> TTTT,KK,r]r], rjrf , K ^irK, a -> KK , r\ir. 

The constraints for decay couplings together with the placement of /o are to linear (n, M 2 ) trajectories rigidly determine 
two nonets with n = 1 and n = 2. The extra state, which was found and investigated in [5,6,8], namely, /Q are (1580), 
should be identified as a scalar glueball, for its decay couplings satisfy the requirements inherent in the gluonium, see 
discussion in [3]. So, 

,/ o bare (1580 ± 40) -» glueball . 
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FIG. 11. The /o-levels in the potential well depending on the onset of the decay channels: bare states (a) 
and real resonances (b). 
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FIG. 12. Linear trajectories on the (n, M 2 ) plane for bare states (a) and resonances (b). 



3. Where is the lightest pseudoscalar glueball? 

In the pseudoscalar sector, in the mass region under discussion, one may expect the lightest pseudoscalar glueball, 
though the opinions about the mass of pseudoscalar glueball are rather different. According to lattice gluodynamic 
calculations, the mass of the lightest pseudoscalar glueball coincides with that of the tensor glueball, that is, it must 
be in the range 2100-2600 MeV [18], while, according to [19], its mass is close to that of the lightest scalar glueball: 
1300-1700 MeV. The plausible existence of the light ^ glueball looks nice, for it might explain a considerable 
production of the ^ states in the radiative J/tp decay, in particular, J/tp — > 777' (according to [20], the admixture of 
the gluonium component in rj' may be rather large, about 10%-20%). However, among narrow resonances one cannot 
see the candidates for pseudoscalar glueball: 77(1295) and ?7(1440) lay on linear qq trajectories (Fig. 4c), though in 
[21] it was estimated that the value of gg component in the 7/(1440) can be not small, 20% ± 20%. Still, it is possible 
that pseudoscalar glueball, after the onset of decay channels, turned into the broad state in the region 1400-1500 
MeV (as it occurred with scalar one) — the experimental data do not contradict this suggestion, see [22] and [23]. 

I am indebted to A.V. Anisovich, D.V. Bugg, L.G. Dakhno, V.A. Nikonov, A.V. Sarantsev for the interest to the 
problem under discussion. The paper is supported by RFFI grant N 01-02-17861. 
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